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The FormCalc package automates the computation of FeynArts amplitudes up to one loop including the 
generation of a Fortran code for the numerical evaluation of the squared matrix element. Major new or 
enhanced features in Version 5 are: iterative build-up of essentially arbitrary phase-spaces including cuts, 
convolution with density functions, and uniform treatment of kinematical variables. The LoopTools library 
supplies the one- loop integrals necessary for evaluating the squared matrix element. Its most significant 
extensions in Version 2.2 are the five-point family of integrals, and complex and alternate versions. 



1. Introduction 



• Enhanced cache functionality, 



FormCalc £Q is a Mathematica package for 
the calculation of Feynman amplitudes. Am- 
plitudes generated by FeynArts [2] are sim- 
plified analytically and converted to a self- 
contained Fortran code for the computation 
of the squared matrix element. The present 
article describes the following features added 
or enhanced in Version 5: 

• Iterative build-up of essentially arbi- 
trary phase-spaces with various cuts, 

• Convolution with density functions for 
external particles, 

• Uniform treatment of kinematic vari- 
ables, 

• Better modularity and code reusabil- 
ity, i.e. less cross-talk between program 
modules. 

The LoopTools library supplies the implemen- 
tations of the scalar and tensor one-loop inte- 
grals necessary for running the Fortran code 
generated by FormCalc. The extensions in 
Version 2.2 are: 

• Addition of the five-point family of func- 
tions, 

• A collective two-point function Bget 
similar to Cget and Dget, 

• Complex versions of most integrals, 

• Run-time selection of alternate versions 
for checking, 

• A command-line interface, 



• Better internal accuracy, and a complete 
quadruple precision build available on 
Intel x86/ifort platforms. 

2. FormCalc 

2.1. Phase-space 

FormCalc 5 no longer uses hand-tailored 
phase-space parameterizations like its prede- 
cessor versions, but builds up the n-particle 
phase-space iteratively by nested integrations 
over the invariant mass Mj and solid angle Qj 
of each outgoing particle i. This procedure is 
encoded in the subroutine Split: 
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Counting the degrees of freedom, there are 
(n — 1) M -integrations and n ^-integrations. 
The corresponding phase-space parameteriza- 
tion is 



1 



•y/s-mi t, 

dMi dfii — 

2y 8 J m2 -\ \-m n 2 

~Mi-m 2 j. 



dM 2 dfl 2 



m^-\ \-m n 



2Thomas Hahn, Michael Rauch 



MPP-2006-7 



x • • • 

/■M n _2-m„_i i 

x / dM„_! dO„_i -^i 

X | dfi„y (1) 

where dfi^ = dcos^ d^,. The particle's mo- 
mentum ki and cos 9i are given in the respec- 
tive decay's rest frame. The (^-integration 
is trivial because of axial symmetry. From 
the practical point of view this looks as fol- 
lows (this code is taken almost verbatim from 
FormCalc's 2to3.F): 

P = 
ex = 
ey = 
ez = 1 

minv = sqrtS 

msum = MASS3 + MASS4 + MASS5 

call Split (5, dble(MASS5), 

& p, ex,ey,ez, minv, msum, fac, 0, 

& Var (XMREM5) , Var (XC0STH5) , Var (TRIVIAL)) 

call Split (4, dble(MASS4), 

k p, ex,ey,ez, minv, msum, fac, 0, 

feVar(FIXED), Var(XC0STH4) , Var(XPHI4)) 

call VecSet(3, dble(MASS3), p, ex,ey,ez) 

One starts with the initial reference direction 
in (ex, ey, ez) and no boost, p = 0. The avail- 
able energy is given in minv and the sum of ex- 
ternal masses in msum. The Split subroutine 
is then called (n — 1) times for an n-particle 
final state. The reference direction, the boost, 
minv, and msum are automatically adjusted 
along the way for the respective remaining 
subsystem and ultimately determine the re- 
maining n-th vector unambigously, which is 
then simply set by VecSet. 

About the integration variables more will 
be said in the next section. For the moment, 
note that the X in XMREM5 refers to the ratio, 
i.e. XMREM5 runs from to 1. The actual inte- 
gration borders are determined internally by 
Split. 

2.2. Variables 

FormCalc 5 introduces a new homogeneous 
system for all (potential) integration vari- 
ables. Each variable is referred to by a pre- 
processor constant, e.g. SQRTS or XCDSTH. The 
following parts can be accessed via preproces- 
sor macros: 



• Var(v) = the actual value of v. 

• Show(v) = the value printed in the out- 
put - to print e.g. t instead of cos#. 

• Lower (w), Upper (v) , Step(w) = the 
lower limit, upper limit, and step width 
of v. If the step is zero, the cross-section 
is integrated over v. 

• CutMin(w), CutMax(u) = the lower and 
upper cuts on v. 

There are two special variables: FIXED for 
fixed values, i.e. no integration, and TRIVIAL 
for trivial integrations. 

2.3. Cuts 

Split allows to place cuts on each M- and 
cos 0- integration. The ip- integration is not 
modified in the present setup. 
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In practice, the application of cuts works as 
follows: 

key = 

CutMin(XMREM5) = E5MIN 

key = key + Cut (CUT_MREM_E, CUT_MIN) 

CutMin(XC0STH5) = -(1 - C0STH5CUT) 
CutMax(XCQSTH5) = +(1 - C0STH5CUT) 
key = key + 

& Cut(CUT_C0STH, CUT_MIN + CUT_MAX) 

call Split (5, dble(MASS5), 

& p, ex,ey,ez, minv, msum, fac, key, 

& Var(XMREM5) , Var(XC0STH5) , Var (TRIVIAL)) 

The value of the cut is deposited in CutMin 
or CutMax and the cut registered by adding 
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an identifier for the cut to the integer key, 
e.g. Cut(CUT_MREM_E, CUT_MIN) specifies a 
cut on the energy (CUT_MREM_£) from be- 
low (C\JT_MIN) which is used to restrict the 
invariant-mass integration (CUT _MREM _E) . 

Note that these cuts really restrict the in- 
tegration limits. They do not introduce veto 
functions (1 in wanted, in unwanted areas) 
into the integrand, which can severely hamper 
convergence. 

In addition, for each external particle i sev- 
eral kinematical quantities are available: 

• momspec(SPEC_M,0 

- mass m, 

• momspec(SPEC_E,i) 

- energy E, 

• momspec(SPEC_K,i) 

- momentum k, 

• momspec(SPEC_ET,i) 

- transverse energy Et, 

• momspec(SPEC_KT,i) 

- transverse momentum fcy, 

• momspec(SPEC_RAP,0 

- rapidity y, 

• momspec(SPEC_PRAP,i) 

- pseudo-rapidity n, 

• momspe c ( SPEC_DELTAK , i ) 

- the difference E — k. 

2.4. Convolution 

With the new variable system, the convolu- 
tion with arbitrary parton distribution func- 
tions can easily be achieved. Three modules 
are already included in FormCalc 5: 

• lumi_parton.F = initial-state partons, 
no convolution. 

• lumi_hadron.F = initial-state hadrons, 
convolution with hadronic PDFs from 
the LHAPDF library |3j. 

• lumi_photon.F = initial-state photons, 
convolution with CompAZ spectrum 0] . 

2.5. Modularity and Code-Reusability 

The choice of parameters is directed as in 
previous versions by the two files process. h 
and run.F, which include one each of 



• Kinematics definitions: 
lto2.F, 

2to2.F, 
2to3.F, 

• Convolution with PDFs: 
lumi_parton . F, 
lumi_hadron . F, 
lumi_photon . F, 

• Model initialization: 
model_sm.F, 
model_mssm.F, 
model_thdm.F. 

There is almost no cross-talk between differ- 
ent modules which are in that sense 'univer- 
sal.' Also, the main program has been rad- 
ically slimmed down in FormCalc 5. It now 
only scans the command line and invokes 

call Processlni ( . . . ) 
call ParameterScan( . . . ) 

All further action is decoupled from the main 
program and can easily be called from any 
application. It is thus straightforward to use 
FormCalc- generated code in own programs. 

3. LoopTools 

LoopTools is a library for the one-loop inte- 
grals. It is based on FF [3] and has a Fortran, 
C/C++, and Mathematica interface. It is ref- 
erenced by the FormCalc-generated code but 
can of course be used also without FormCalc. 

3.1. Five-point functions 

The most significant addition is the five- 
point family of functions - the scalar inte- 
gral Eg and the tensor coefficients up to rank 
four. The default versions use the Denner- 
Dittmaier decomposition in which one in- 
verse power of the Gram determinant is can- 
celled and thus has considerably better nu- 
merical stability than the Passarino-Veltman 
|7] decomposition, which is also available for 
comparison. 

3.2. Combined two-point function 

The two-point functions have been united 
into the Bget function which works similar to 
its Cget, Dget, and Eget counterparts, in par- 
ticular it caches its results. Compatibility rou- 
tines for the old BO, Bl, etc. are of course avail- 
able. The reason is mainly cache efficiency in 
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view of the five-point decomposition: 
Eget ^5 Dget 

^(5- 4) Cget 

(5 • 4 • 3) Bget (2) 

3.3. Complex versions 

Versions of the LoopTools functions for 
complex parameters have been added as far 
as they are contained in FF, i.e. currently 
only some special cases for the complex Do 
are available. 



real versions 


complex versions 


AO, A00 
BO, Bl, ... 
BOi, Bget 
CO, COi, Cget 
DO, DOi, Dget 
EO, EOi, Eget 


AOC, AOOC 
BOC, B1C, . . . 
BOiC, BgetC 
COC, COiC, CgetC 
DOC, DOiC, DgetC 
EOC, EOiC, EgetC 



3.4. Alternate versions 

For some functions alternate versions exist, 
most of which are based on an implementation 
by Denner. The user can choose at run-time 
whether the default version 'a' (mostly FF) 
or the alternate version 'b' (mostly Denner) 
is used and whether checking is performed. 
This is determined by the version key: 

0*key compute version 'a', 

l*key compute version 'b', 

2*key compute both, compare, return 'a', 

3*key compute both, compare, return 'b'. 

Usage is as in 

call setversionkey(2*KeyC0 + 3*KeyD0) 

The following keys for alternate versions are 
currently available: KeyAO, KeyBget, KeyCO, 
KeyDO, KeyEget, KeyEgetC. KeyAll comprises 
all of these. 

3.5. Debugging 

Debugging output can be turned on simi- 
larly with e.g. 

call setdebugkey (DebugC + DebugD) 

Identifiers range from DebugB to DebugE and 
are summarized by DebugAll. The integrals 
are listed in the output with a unique serial 
number. If the list of integrals becomes too 
long, one can select only a range of serial num- 
bers for viewing, as in 

call setdebugrange(4711, 4715) 

This makes it easy to monitor 'suspicious' in- 
tegrals. 



3.6. Environment variables 

Version key, debug key, and debug range 
can be set also from the 'outside,' i.e. through 
environment variables, thus making it unnec- 
essary to recompile the code. Unfortunately, 
the Key J and Debug!" labels cannot be used 
on the shell command- line, but -1 offers a con- 
venient way to set all bits of an integer, thus 
having the same effect as KeyAll or DebugAll. 
For example: 

setenv LTVERSION -1 
setenv LTDEBUG -1 
setenv LTRANGE 4711-4715 

3.7. Command-line Interface 

The new command-line interface is useful 
in particular for testing and debugging. It 
lists the iV-point scalar and tensor coefficients 
corresponding to the number of arguments, 
i.e. 3 arguments — B, 6 arguments = C, 
etc. Currently only the real versions are ac- 
cessible through the command-line interface, 
mainly because it is not straightforward syn- 
tactically to specify complex parameters on 
the command-line. 

3.8. LoopTools caches 

The internal caching mechanism is mean- 
while used by the following functions: Bget, 
Cget, Dget, Eget, BgetC, CgetC, DgetC, and 
EgetC. Obviously the former system with 
its getcachelast and setcachelast call for 
each individual cache was no longer practica- 
ble for all those caches. 

The new cache-management functions op- 
erate on all caches simultaneously: 

• call clearcache 

- clears all caches, 

• call markcache 

- marks the current position, 

• call restorecache 

- reverts to the last marked position. 

For compatibility, LoopTools still includes two 
routines getcachelast and setcachelast. 
They work only approximately as before, how- 
ever, and are therefore deprecated. 

Furthermore, the cache mechanism itself 
has been improved and now uses a binary 
search method. Cache lookups are thus now 
faster, but as this has not often been a bot- 
tleneck, the impact on program performance 
will typically be minor in most applications. 
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3.9. Accuracy 

The accuracy of the tensor reduction has 
been improved through the use of an LU de- 
composition for the Gram matrix and quadru- 
ple precision at strategic points internally. 
Thanks to the Intel ifort compiler, the latter 
is now widely available on the x86 platform. 

A complete quadruple-precision build can 
be chosen with an alternate makefile. In par- 
ticular the quadruple-precision version of the 
frontend can be useful to check cases where a 
severe loss of precision is suspected. 

4. Summary 

FormCalc 5 is the current release of the 
FormCalc package with the following new fea- 
tures and extensions: 

• New kinematics routines can build up 
almost arbitrary phase-spaces with a 
wide range of cuts possible. 

• FormCalc's driver code consists of in- 
tuitively programmed, concise modules 
with almost no cross-talk between them. 
Calling FormCalc-generated code from 
other applications has been much sim- 
plified. 

• Convolution with arbitrary distribution 
functions is possible. The following 
common cases are available out-of-the- 
box (if FormCalc came in one): 

— partons (no convolution), 

— hadrons (uses LHAPDF), 

— photons (uses CompAZ). 

• All possible integration variables are 
treated uniformly and detailed kinemat- 
ical information is available for all exter- 
nal legs. 

LoopTools 2.2 is the latest version of the 
LoopTools library with the following novel- 
ties: 

• The scalar five-point function including 
its tensor coefficients up to rank four is 
provided. 

• Complex versions are available as far as 
implemented in FF. 



alternate versions, by program call or 
through an environment variable. 

• In the same manner, debugging output 
can be turned on and off individually at 
run-time. 

• The cache system has been extended 
and improved. 

• A new command-line interface is useful 
for testing and debugging. 

• The accuracy of the tensor reduction has 
been enhanced. 

• A complete quadruple-precision build 
can be chosen with an alternate makefile 
(e.g. make -f makefile . quad-if ort). 

Both packages are available as open source 
and stand under the GNU Lesser General 
Public License (LGPL). They can be obtained 
from the Web sites 

http : / / www . f eynarts . de/f ormcalc, 
http : / / www . f eynarts . de/looptools. 
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• Checking can be enabled at run-time, 
individually for all integrals which have 



